Search for the second forbidden beta decay of to the ground state of ^Be 
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A significant decay branch of *B to the ground state of *Be would extend the solar neutrino spec- 
trum to higher energies than anticipated in the standard solar models. These high-energy neutrinos 
would affect current neutrino oscillation results and also would be a background to measurements 
of the hep process. We have measured the delayed alpha particles from the decay of *B, with the 
goal of observing the two 46-keV alpha particles arising from the ground-state decay. The *B was 
produced using an in-flight radioactive beam technique. It was implanted in a silicon PIN-diode 
detector that was capable of identifying the alpha-particles from the *Be ground state. From this 
measurement we find an upper limit (at 90% confidence level) of 7.3 x 10~^ for the branching ratio 
to the ground state. In addition to describing this measurement, we present a theoretical calculation 
for this branching ratio. 

PACS numbers: 23.40.-s, 23.40.Bw, 23.60.-|-e, 26.65.-|-t, 27.20.+n 



I. INTRODUCTION 

There have been several experiments [H, 0; 13 done re- 
cently to measure the alpha-particle spectrum following 
the allowed (2+ to 2+) beta decay of ^B to the broad 
first excited state at 3 MeV in ®Be. The neutrinos from 
this decay provide the signal for solar neutrino measure- 
ments, such as SNO and SuperKamiokande, and precise 
knowledge of their spectrum at creation, combined with 
accurate measurements in the solar neutrino detectors, is 
important for understanding neutrino oscillations. The 
only higher energy neutrinos from the Sun come from the 
hep [1] process (weak capture of a proton by '^Hc produc- 
ing e'^ + v), expected to be at a rate roughly 10^'^ of the 
®B neutrinos, and from the second forbidden (2+ to 0^) 
decay of ®B to the ground state of ^Be. Specifically these 
reactions are 



'He 



^''He-^e^ 



and 



^B 



*Bc 



18.77 MeV 



16.95 MeV. 



An ft value for the latter decay, estimated using 
log^o ft = 13.3 obtained from the 2+ to 0+ decay of ^^Cl, 
suggests the *B 2+ to 0+ decay should have a branching 
ratio of order only 10"^. The ^B decay has an unusually 
large energy, and second-forbidden decays should have 
an extra factor of /c^ (k is the momentum transfer to the 
leptons) in their decay rate. Comparing the square of 
the endpoint energy for ^B with that for '^^Cl suggests 
an additional factor of 570 in the rate for the second 
forbidden decay of ^B. This rate would correspond to 



•Presently at Nuclear Safeguards Science and Technology, Los 
Alamos National Laboratory, Los Alamos, NM 87544 
t Presently at Department of Physics and Astronomy, University of 
British Columbia. Vancouver, EC V6T IZl 



a log;^o of 10.5, a value still in the expected range for 
second-forbidden transitions. The corresponding branch- 
ing ratio would be 3 x 10~^. If the branching ratio were 
as large as about 1%, it would contribute significantly to 
the neutrino spectrum observed at SNO. 

Decay to the ground state of ®Be is difficult to mea- 
sure, since that state decays to two alpha particles pro- 
ducing only 92 keV. We are unaware of previous attempts 
to measure this branch. Published estimates range from 
3x 10-6 d to 0(10-4) The latter estimate includes 

only the vector contribution to the second-forbidden de- 
cay. In this paper we describe the technique of our mea- 
surement, give the result, and present theoretical calcula- 
tions based on the shell model and on measured gamma- 
ray transition matrix elements. 



II. EXPERIMENT 

We made a ^B beam using an in-fiight production tech- 
nique. The ®B were implanted in a silicon detector that 
was capable of measuring the energy from the delayed al- 
pha decay of the ^Be ground state, as well as from most 
of the spectrum of delayed alpha decay of the ^Be first 
excited state. The decays were tagged by detecting the 
decay positrons in a scintillator. The silicon detector was 
calibrated using back scattered a particles of known en- 
ergy and also using delayed alpha decays of implanted 



A. Radioactive beams 

The *B for implantation in the detector was made us- 
ing the ■^He(6Li,8B)n reaction. A 32-mm-diamcter gas 
cell was located at the focal point of the analyzing mag- 
net of the University of Washington Tandem accelerator. 
This gas cell had entrance and exit windows of 2.5-/zm 
Havar H and contained ^He at about 0.6 bar. The inci- 
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dent beam was 24-MeV ^Li, which, after accounting for 
energy loss in the windows and gas, produced *B in the 
forward direction with an energy centered at 15.5 MeV 
and a spread of 1.1 MeV, resulting from energy loss in 
the gas. 

The beam transport following the analyzing mag- 
net consisted of a 30° switching magnet followed by a 
quadrupole doublet. These elements are normally used 
to focus the analyzed beam from the accelerator onto the 
target in the scattering chamber. For the present exper- 
iment we had to optimize the capture and transport of 
the radioactive beam through the 20-mm aperture for 
our detector, which was located at the target position in 
the scattering chamber. We found, using transport 
[§| and TURTLE that this could be done by first 

tuning the beamline following the gas cell appropriately 
for transport of a straight-going monoenergetic beam of 
the desired magnetic rigidity, and then by increasing the 
quadrupole strength by about 10%. The energy accep- 
tance of this system was only about 2%, which is less 
than the energy spread resulting from the gas cell, so 
there was no reason to have either a bigger cell or denser 
gas. We tuned the radioactive beam transport elements 
by tuning first a 7.44-MeV ^Li'^^ beam (matching the 
15.5-MeV *B^+ radioactive beam in magnetic rigidity) 
through a small aperture. Then, while irradiating the 
gas cell with the 24-MeV ^Li beam, we would monitor 
the ®B implantation rate and increase the quadrupole 
current to maximize the rate. We also swept the switch- 
ing magnet field, which we monitored with a Hall probe, 
to make sure we were accepting the central energy from 
the gas cell. 

As the rigidity of the primary beam is higher than that 
of the radioactive beam, degraded primary beam is a seri- 
ous contaminant. In order to minimize degradation of the 
primary beam after the analyzing magnet, we regulated 
the accelerator voltage with the Generating Volt Meter 
Regulator [T]| and completely retracted the beam regu- 
lation slits located downstream of the analyzing magnet. 
Particular attention had to be paid to the beam tuning, 
in order to minimize the amount of beam striking other 
limiting apertures in the beamline. 

In addition to the *B beam, we produced a ^^N beam 
for calibration of the silicon detector. This beam was 
made using the 3He(i°B,i2N)n reaction, with a 35-MeV 
^"^B beam, producing 24.2-MeV ^^N^+, in an analogous 
manner. The ^^N energy was chosen to produce the same 
implantation depth for the ^^N as for the ^B. The choice 
of the ^■^N charge state was determined by the back- 
ground, because the most-probable charge state of ^^N, 
6-I-, was too close in rigidity to the ^"^B primary beam. 
Using the coincidence with the scintillator, we were able 
to tag about 3 — 5 ^B per second with 170 to 330 particle 
nA of "^Li, and about 0.5-1 ^^N per second with 50 to 100 
particle nA of ^°B. For the *B beam, the background in 
the silicon detector, primarily consisting of degraded ^Li 
from the primary beam, was 0.6 — 1 kHz. For the ^^N 
beam, the degraded ^°B background was 2 kHz. This 



background rate in our detector set an intensity limit for 
the primary beam. 



B. Detector 

The detector assembly consisted of a 500- /im thick by 
18-mm square Hamamatsu S-3204-06 PIN diode [12.], in 
front of which was a 51-mm cube of Bicron-400 scintil- 
lator with a 20-mm-diameter hole on the beam axis. A 
lead shield 6. 4- mm thick with a 17-mm-diameter aperture 
was placed in front of the scintillator, and the scintillator 
was coupled on the side to a photomultiplicr tube. This 
thickness of lead was sufficient to absorb positrons from 
decaying ^B that stopped on it. Ions that passed through 
the lead aperture would stop in the PIN diode, and were 
not able to hit the sides of the hole in the scintillator. 

The 15-MeV ^B ions were implanted approximately 20 
/im into the silicon, and about 20% of the decays pro- 
duced positrons that would enter the scintillator, tag- 
ging the decay. In order to minimize distortions in the 
a-particle spectrum from positron energy deposition, the 
scintillator detector was placed in front of the silicon de- 
tector, thus minimizing the track length of the tagged 
positrons in the silicon detector. Tagging with positrons 
that passed through the silicon detector was not practi- 
cal, as their energy deposition would be larger than the 
energy from the decay from the ^B ground state. Indeed, 
as is discussed in section III D 21 positrons that scatter 
through large angles affect the efficiency of the detec- 
tion. The PIN diode was mounted on a custom charge- 
sensitive preamplifier Ts*! built on a 5-cm-square circuit 
board. The preamp and detector were cooled by a ther- 
moelectric cooler coupled to a 19-mm-wide and 3-mm- 
thick copper strip that ran the height of the preamplifier 
circuit board. This assembly is shown in Fig. [1] 

When cooled to 10° C, the measured leakage current 
of a new PIN diode was typically 10—12 nA. The resolu- 
tion (FWHM) was 7.8 keV for 81-keV gamma rays, and 
a threshold about 30 keV was practical. During experi- 
mental runs, the reverse current would rise to 30 — 60 nA 
after 2 weeks of irradiation, at which point the resolution 
would decrease to approximately 10 keV at 81 keV. 



C. Electronics and data acquisition 

The photomultiplier on the scintillation detector pro- 
duced both a timing and a pulse-height signal. The tim- 
ing signal started a time to amplitude converter (TAG). 
The output of the charge-sensitive preamplifier for the 
silicon diode detector was split into two channels, each 
with an Ortec 572 amplifier. The maximum range of the 
high gain channel corresponded to about 1 MeV, while 
the maximum observable energy in the low gain chan- 
nel was 8 MeV. The two unipolar signals were digitized, 
as was the pulse-height signal from the scintillator. The 
bipolar outputs of the amplifiers for the silicon detec- 
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Also, the irradiations were stopped at hourly intervals, 
and a measurement of the pulser spectrum for vari- 
ous settings of the pulser attenuation circuit was made 
(1.0x,0.5x,0.2x,0.1x). These measurements allowed 
for the monitoring of both gain and linearity as a function 
of irradiation time. 

From the event-mode data, singles histograms of the 
silicon-detector energy spectra, the TAG signals, and the 
scintillator pulse height were formed. Prompt coincident 
events were retained by applying cuts (discussed in sec- 
tion |lll| to the TAG and scintillator spectra. 
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D. Calibration 



1. Silicon detector energy calibration 



FIG. 1: Schematic of the detector assembly. The various hole 
diameters are given in the text. The cooled copper mount 
is not shown, but was located behind the preamplifier circuit 
board. 



tor generated timing signals, which were combined in an 
OR module, and the resulting logic signal was delayed 
and used to trigger the data acquisition system (DAQ), 
contained in a single computer-automated measurement 
and control (GAMAG) crate. A parallel output of the 
OR unit was used to stop the TAG, which was started 
by the signal from the scintillator. This TAG signal was 
also digitized. Thus any signal over the 30-keV threshold 
from the silicon detector would trigger the DAQ, while a 
coincidence (within a 5-/xs time window) would also pro- 
duce a TAG output. The threshold on the scintillator 
was set high enough that its rate was substantially lower 
than that in the silicon diode detector. The DAQ system 
was operated in event mode, with two silicon detector en- 
ergy signals, one scintillator pulse height, and one timing 
signal per event. 

A precision pulser, coupled to the silicon detector 
preamplifier, was used for monitoring dead time and am- 
plifier gain shifts. A scaler module was used to count 
pulser triggers, pulses from the beam-current integrator, 
the number of TAG outputs, and the number of DAQ 
triggers. The pulser and the electronics feeding the DAQ 
system were located in a temperature-stabilized rack, 
which was controlled to within ±0.1°G. During regular 
data taking the pulser was set to correspond to about a 
0.5 MeV signal in the silicon detector. The DAQ GA- 
MAG crate was located in the temperature-controlled 
electronics rack and was controlled by a Wiener GG32 
peripheral component interconnect (PGI) crate controller 
[l4| . For data acquisition, initially we used the JAVA- 
based jam [l^ system, but subsequently switched to the 
Objective-G-based ORCA system. 

The gain stability of the silicon-detector preamplifier 
and associated electronics was monitored by observing 
the pulser signal in the high- and low-gain channels. 



The energy calibration for the ®B(2+) **Be(0+) de- 
cay measurement was obtained from a ^'^^Ba radioactive 
source, which produces a number of low-energy gamma- 
ray lines between 80 — 400 keV. These photons convert by 
Gompton scattering and the photoelectric effect in the sil- 
icon. In contrast, measurement of the *B(2+) — > ®Be(0"'") 
decay rate involved the detection of two low-energy alpha 
particles, which directly ionize the detector medium. The 
ionization energy for gamma rays in silicon, 3.68 ± 0.02 
eV, is larger than the ionization energy for alpha par- 
ticles in silicon, 3.62 ± 0.02 eV. This difference, a ra- 
tio of 1.017 ± 0.008, must be accounted for when relat- 
ing a calibration obtained from gamma-rays to the mea- 
sured alpha-particle energy spectrum [iTj . In addition, 
the fractional loss of energy to non-ionizing processes in- 
creases as the alpha particle incident energy decreases; 
this pulse-height deject goes undetected by the silicon de- 
tector. 

Tabulated values of the expected nuclear and total 
stopping powers have been published in [l^. At low en- 
ergies, these tabulated nuclear stopping powers are calcu- 
lated via a model, as available experimental data are lim- 
ited. In order to determine the expected pulse height of 
the ®B(2+) ®Be(0"'") decay signal reliably, it was neces- 
sary to measure the pulse-height defect of alpha particles 
in our silicon detector. 

The uncertainty in the ionization energy ratio affects 
the determination of the pulse-height defect, but, as will 
be seen below, the calibration of the low-energy alpha 
particles is not affected by it. Thus we divide the cal- 
ibration obtained from ^'^^Ba by 1.017 to get the cor- 
responding calibration for alpha particles without the 
pulse-height deficit, and these energies are used for our 
calibrations. The pulse-height deficit (with uncertainties) 
is applied subsequently. 

In order to measure the pulse-height defect of low- 
energy alpha particles in silicon, '^He ions were produced 
at eight different energies from 90—475 keV using the UW 
Tandem accelerator with the Terminal Ion Source. These 
alpha particles were backscattered from a gold monolayer 
evaporated on a 50-/xg/cm^ carbon target foil into the 
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collimated silicon detector located at 100 degrees. As 
the energy the particles deposit in the active part of the 
detector is reduced from the incident energy by the Si02 
deadlayer on the front of the detector, it was necessary 
to quantify the energy loss in the deadlayer. This was 
done by repeating the measurements for each incident 
energy with the detector face angled at numerous posi- 
tions between —45 and +45 degrees to the scattering axis. 
At each angle, the effective thickness of the deadlayer 
through which the particles passed is Ace/ cos 0, where 
Aa; is the physical deadlayer thickness. The value of Aa; 
was determined by combining all of the measured results 
and performing a global fit to the measured energy (as 
determined from the ^'^^Ba calibration, adjusted by the 
ratio of ionization energy for electrons and alpha parti- 
cles) vs. detector angular position. From this series of 
experiments, the pulse-height defect was found to vary 
from 7.9±0.3 to 11.4±0.7 keV for alpha-particle incident 
energies in the active part of the detector varying from 
55.9 to 391.4 keV. The measured energy loss in the dead- 
layer corresponded to 43.1 ± 0.4/ig/cm^ of Si02. There 
is an additional systematic uncertainty of 0.8% in these 
pulse-height defects, arising from the uncertainty in the 
ratio of ionization energy for alpha particles to that for 
photons [13] ■ 

The low-energy alpha calibration, including corrections 
for pulse-height defect and for the ratio of ionization en- 
ergy for electrons and alpha particles, was cross-checked 
using the decay of implanted ^^N, which beta-decays to 
the 7.65-MeV state of ^-^C, which then decays to three 
alpha particles with a total energy release of 379 keV. 

Two of the alpha particles come from the breakup of 
the ground state of *Be. By comparing the known energy 
of this decay (accounting for decay kinematics, discussed 
in the following section) to the detected energy deposi- 
tion, we were able to determine the pulse-height defect 
independently, assuming that the energy dependence of 
the pulse-height defect was that given by the predicted 
non-ionizing stopping power . The spectrum obtained 
from the ^^N decay is shown in Fig. [2l The results of the 
^^N decay and the single alpha-particle measurements are 
shown in Fig. [31 along with the predicted pulse- height de- 
fect obtained from [18]. The pulse-height defect obtained 
from the single alpha-particle measurement is calculated 
to be 7.5 ± 0.3 keV at 46 keV. On the other hand, the 
pulse-height defect from the ^^N decay measurement is 
found to be 4.5 ± 0.9 keV at 46 keV. We use these two 
values as lower and upper bounds when constructing the 
energy window for the expected signature from *B decay. 

As the low-energy alpha calibration includes correc- 
tions for both pulse-height defect and for the ratio of 
ionization energies, and as the former is obtained using 
the latter, uncertainties in this ratio will cancel in the 
calibration for the *B measurement. The photons are 
just an intermediate calibration. The determination of 
the pulse-height defect links a known alpha-particle en- 
ergy to a photon energy, and then the use of that photon 
for calibrating the measurement of *B establishes that 
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FIG. 2: The ^'^N 3a decay peak observed in the high gain 
channel, fit with a Gaussian peak on a linear background. 
The energy deposited, including recoil and positron summing, 
is predicted to be 388.5 ±1.2 keV from our Monte Carlo sim- 
ulations. The center of the peak appears at 372.5 ± 3.0 keV, 
indicating a total 16 ± 3 keV pulse-height defect. 
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FIG. 3: Pulse-height defect (PHD) measured using single al- 
pha particles (solid circles) and from ^^N decay (open circles). 
The shaded band is the PHD value obtained from integrating 
the stopping-power tables in p^ . including a 10% uncertainty. 



calibration in terms the known alpha-particle energy. 

The low-gain channel, which contained the 
^B(2-|-) ^^Be(2-f) decays, was calibrated using the two 
lines from ^^N decay, resulting from the 3a decay of the 
levels at 7.65 and 12.71 MeV. These are illustrated 
in Fig. m The energies of these lines were adjusted for 
positron summing, ^^C recoil and pulse-height defect. 



2. Recoil and energy deposition calculations 

Two corrections are required for the analysis of both 
^B and ^^N energy spectra: the recoil energy imparted 
by the leptons to the daughter nucleus, and the amount 
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FIG. 4: The N spectrum observed in the low-gain chan- 
nel. The peaks near channels 80 and 1300 correspond to total 
energies from 3a decay of 388. 5± 1.2 and 5443±6 keV, respec- 
tively, including recoil and positron summing. The continuum 
is dominated by positrons backscattering from the detector 
mount and thus passing through the detector twice. In ^^N 
decay, only 3% of the positrons are accompanied by alphas. 



of energy deposited in the silicon detector by the out- 
going positron. As a result of the high decay energy, a 
significant amount of recoil energy (up to 19 keV) can 
be imparted to the ^Be daughter nuclei by the outgoing 
leptons, resulting in an overall broadening of the energy 
window in which the decay signal is expected to lie. The 
^^C daughters receive less recoil energy (up to 3.3 keV). 
Using the predictions of the positron-neutrino correla- 
tion from the calculations described in Sec. IIVI we found 
a distribution of recoil energy between and 19 keV for 
the ^B(2+) *Be(0"'') decay, which was peaked near the 
maximum energy. The Gamow- Teller decay of "'^^N yields 
a fairly uniform distribution of recoil energy between 
and 3.3 keV. 

The additional energy deposited by the escaping 
positron adds to the decay energy and results in an over- 
all energy shift. We minimized the magnitude of this 
energy shift by implanting the *B and ^^N ions near the 
upstream surface of the silicon detector. The distribu- 
tion of the positron energy deposition was computed by 
Monte-Carlo simulation, using measured stopping powers 
from [Hi. The ®B and ^^N implantation depths, 20.2±0.3 
and 18.7 ± 1.5 microns respectively, were obtained from 
p^ . The simulated energy distributions were fit with a 
Landau distribution, and the most likely energy deposi- 
tion for positrons coming from ^B(2+) ^Be(0'^) and 
-^^N decay was 9.5 ± 0.9 and 7.8 ± 0.6 keV respectively. 

A significant number of positrons undergo large an- 
gle scattering, in either the silicon or the copper detec- 
tor mount, before reaching the scintillator and result- 
ing in a coincidence trigger. Their energy deposition 
in the silicon detector is large enough to remove the 
event entirely from the expected energy range for the 
®B(2+) — > ^Be(0+) decay. As described below, we deter- 



mined that positron summing reduced our efficiency for 
detecting ground state decays by 11.5%. 



III. DATA ANALYSIS 

Following a series of commissioning runs, we accumu- 
lated 2.5 X 10^ ^B decay events over the course of two 
experimental runs. After selecting prompt events and 
applying cuts on the amplitude of the pulse in the scin- 
tillation counter, 2.0 x 10^ events remained. The follow- 
ing discussion applies to both runs and is illustrated by 
examples from the second. 

Following each run, the data in each of the silicon de- 
tector channels were corrected for gain drifts (the aver- 
age correction was less than 0.3%), and the high-gain 
channel was calibrated using the information from the 
^^•^Ba calibration data. A correction was also made for 
the alpha/gamma ionization energy difference. No ad- 
justments or calibrations were made to the scintillator or 
TAG spectra. The events of the final data set consisted of 
the corrected silicon-detector values and the correspond- 
ing raw scintillator and TAG values. This processing was 
applied separately to the data obtained using the ^He 
and *He target gas. The ratio of random to true coinci- 
dence rate, in the energy region where the ground state 
decay is expected, was found to be approximately 0.5% 
for both the "^He and ^He data. Since these data were 
subtracted to obtain the final data set, we made no fur- 
ther correction for accidental coincidences. 

The spectrum from the low-gain channel, which is 
dominated by the *B(2-|~) — >*Be(2-|-) decay, is shown in 
Fig. [51 Above about 8 MeV, Monte Garlo simulations in- 
dicate the spectrum is distorted by alpha particles that 
escape from the 20-/im implantation depth. Additional 
small distortions result from the scintillator threshold, 
which requires positrons over about 1.5 MeV, and conse- 
quently favors lower energy alpha particles. Finally, sig- 
nificant distortions result from backscattered positrons, 
many of which come from the copper cooling mount be- 
hind the detector. Presumably as a result of these distor- 
tions, which are examined in detail in Ref. (20| . we were 
unable to get a useful R-matrix fit to this spectrum. 

To determine the energy window in which the 
^B(2+) — + ^Be(0+) decay should be located, we begin by 
adding the recoil energy, distributed between and 19.3 
keV, to the total kinetic energy (91.8 keV) of the two 
alpha particles arising from ^Be decay. To calculate the 
pulse-height defect of these alpha particles, we assume 
that the total energy is shared equally between the two 
alpha particles, thus giving a minimum and maximum ki- 
netic energy for each alpha particle of 45.9 keV and 55.6 
keV, respectively. The range of the pulse-height defect 
for these alpha particles was measured to be 4.5—7.5 keV. 
Thus, before accounting for the positron energy loss, the 
ground state decay could be expected to appear between 
9 and 15 keV lower than the sum of the decay and recoil 
energies. We convolve the calculated alpha-particle en- 
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FIG. 5: Spectrum from the low-gain channel of the silicon 
detector, dominated by the *B(2-|-) ^*Be(2-|-) decay. As 
discussed in the text, the 20- implantation depth limits 
the reliability of this spectrum to energies below 8 MeV. 



ergy distribution with a Gaussian detector response with 
the typical 9-keV FWHM to produce an expected energy 
spectrum. The pulse height defect shifts this spectrum 
downward between 9 and 15 keV. For each of these cases 
we convolve the shifted spectrum with the spectrum of 
positron energy deposition calculated using the Pene- 
lope Monte-Carlo package [2l[. That spectrum corre- 
sponds to positrons emitted isotropically but depositing 
over 2.5 MeV in the scintillator. The spectrum has a peak 
at about 9 keV, corresponding to positrons that exited 
the silicon in the direction of the scintillator, but there 
is a long tail produced by positrons that either backscat- 
tered from the copper mount or that passed through sig- 
nificant thickness of silicon before a large-angle scatter- 
ing into the scintillator. The final spectra, for the two 
limiting values of pulse height defect, including positron 
summing, are illustrated in Fig. [51 We select the en- 
ergy range from 75 to 120 keV for the region of interest. 
This region contains between 88.1% and 88.8% of the to- 
tal events, depending on the assumed value for the pulse 
height defect. 

In this energy window the rates of the two background 
processes in ^Be, gamma-decay of excited states of ^Be 
and alpha particles from the broad ^B(2+) — > *Be(2+) 
decay spectrum, were calculated to be two orders of mag- 
nitude below the sensitivity of this experiment. 

The spectra obtained in the second experimental run, 
with '^He (signal) and "'He (background) in the gas cell, 
are shown in Fig. [T] The most striking feature is the 
large broad peak of coincident events, which was present 
in our data regardless of the target gas, and (at a slightly 
lower rate) even without the beam. The rate of events 
in this peak, 20 mHz, was relatively constant over a long 
time span (^ 4 years) , regardless of the time between ir- 
radiations in the scattering chamber. We speculate that 
this peak is due to minimum-ionizing events from am- 
bient radioactivity which pass through both the silicon 
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FIG. 6: Monte Carlo simulations of silicon detector response 
to decays to the *Be ground state for coincident events in 
scintillator and Si detectors. The initial alpha-particle energy 
spectrum is modified by recoil and positron energy summing. 
The two curves correspond to the two extreme values for the 
pulse height defect. The arrows indicate the region chosen for 
searching for ground state decays. 
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FIG. 7: Measured signal (solid) and background (dashed) co- 
incidence silicon-detector spectra for the high-gain channel 
with scintillator threshold applied. 



and scintillator detectors. However, increased shielding 
and the addition of scintillator veto counters were not 
effective in reducing the rate. 

Analysis of these events shows that they are dominated 
by low-energy scintillator events, and we can improve 
our experimental sensitivity by optimizing the scintillator 
minimum threshold. The optimum scintillator threshold 
was obtained by finding the minimum uncertainty in the 
difference of the time-normalized signal and background 
counts between 75 — 120 keV as a function of scintilla- 
tor threshold. We maintained the average beam current 
in the measurements with '^He and ^He to be identical 
within 1%, so normalizing by time was equivalent to nor- 
malizing by integrated current. We subtracted the nor- 
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malized background given by the ^He data from the "^He 
data. 

Along with the background discussed above, there was 
additional background present in the spectrum measured 
using the signal (^Li + '^He) reaction which was not 
present in the background (^Li + ''He) measurement. 
We hypothesize that this additional background resulted 
from positrons emitted from *B implanted on surfaces 
other than the silicon detector. These positrons may 
pass through both the silicon and scintillator detectors 
and result in a false coincidence trigger and a continuum 
of events in the energy spectrum. 

Since this background had the same general shape as 
the larger background observed with the beam off or with 
''He in the gas cell, we assume it is composed primarily of 
a spectrum of the same shape, but we will allow for the 
possibility of additional background in the energy range 
of interest. Such a background would be identified by 
additional counts in the spectrum adjacent to the energy 
window for the ground state decay. 

In order to account for the first part of this background, 
we parameterized the shape of the background spectrum 
by fitting the measured background energy spectrum 
from 50—400 keV, as shown in the first plot of Fig. [HI The 
parameterization function consisted of a Landau distribu- 
tion and a linear term. Holding the parameters from the 
background fit fixed and allowing for a free multiplicative 
normalization factor, the parameterization function was 
then fit to the signal energy spectrum over an energy win- 
dow from 140 — 400 keV, as shown in the second plot of 
Fig.m Excluding the area around the ®B(2+) ^Be(0+) 
decay signature window ensures that the normalization 
was not influenced by counts in that window. The num- 
ber of expected background events in the signal region 
was obtained by integrating the fit function over the en- 
ergy region of interest. In order to identify the second 
part of this background, the excess of counts above this 
calculated background was computed in the signal win- 
dow (75 — 120 keV) and two adjacent energy windows 
(50 - 75 keV and 120 - 140 keV). 

The excess of net counts above the fit background in all 
three energy windows confirms the second background, 
which is not present in the ^He measurements from which 
the initial background estimates are derived. To account 
for this additional background, we interpolated between 
the two adjacent windows to determine the number of 
background events in the primary energy window. We 
then subtracted this background from the total number 
of excess events obtained from the difference between the 
data and the fit, to get the net excess of signal events. To 
account for the loss of efficiency due to positron backscat- 
tering, we multiplied the net excess of counts by 1.13. 
The branching ratio was calculated by dividing the cor- 
rected net excess number of counts by the total number 
of detected coincident events. The results of this analysis 
are given in Table [J 

The data from the two experimental runs, given in 
tableHl combined with the total of 1.96 x 10^ decays col- 




Energy in silicon detector (keV) 




Energy in silicon detector (keV) 



FIG. 8: Data from the second experimental run. The top plot 
is the data measured with *He in the gas cell, with the fit to 
this background. The bottom plot is the data measured with 
^He in the gas cell, with a fit to the data above 140 keV using 
the parameters obtained from the background spectrum, but 
allowing for an adjustable multiplicative normalization. The 
dashed range of the second curve denotes the energy region 
of interest for signal extraction. 



lected, yields a branching ratio for the ^B(2+) **Be(0"'") 
decay of (2.6 ± 3.3) x 10~^. This result is consistent with 
no measured signal and corresponds to an upper limit on 
the branching ratio of 7.3 x 10~^ at a 90% confidence 
level. 



IV. THEORETICAL ESTIMATE OF 
BRANCHING RATIO 

One published estimate of the 8B(2+) ^ «Be(0+) 
decay rate is based on conserved vector current (CVC), 
which relates the vector contribution to the beta-decay 
rate to the gamma decay of the 16-MeV isobaric analogue 
state of ^Be [12]. These calculations do not consider the 
contribution of the axial matrix element to the decay 
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TABLE L Number of measured events and predicted background in the spectrum obtained with '^He gas in the cell. Both 
experimental runs are combined. 



energy window (keV) 


events from data 


bkg. events from fit 


net remaining 


interpolated bkg. 


net excess 


50-75 


282 ± 16.8 


255.7 ± 16.0 


26.2 ±23.2 






75-120 


661 ±25.7 


575.4 ± 24.0 


85.6 ±35.2 


40.4 ± 44.4 


45.2 ± 56.7 


120-140 


589 ± 24.2 


573.9 ± 23.8 


15.1 ±34.0 







rate. Another estimate [5] is based on several different 
model calculations of the axial-vector second-forbidden 
contribution to the decay to the first excited state of 
^Be. Those authors suggest that, based on the similarity 
of the intrinsic wave functions of the ground and first 
excited states of *Be, these calculations can be applied 
to the ground state decay, and they find a branching 
ratio of approximately 3 x 10~^. Hovi^ever, since there 
is no experimental information available to constrain the 
contribution of the axial matrix element to the transition, 
a strong axial contribution cannot simply be ruled out. 

In this section, we discuss our a ppli cation of the frame- 
work developed by Walecka [23l. |24|| and Donnelly and 
Haxton [2^ for calculating the rates of electromagnetic 
and semileptonic processes. The details of the cal- 
culations discussed in this section are specific to the 
®B(2+) — > ^Be(0"'") decay. Our application of this frame- 
work is described in detail in Ref. [20| . 

A. Shell model estimate 

The rate of the ^B(2+) -> 8Be(0+) second-forbidden 
^B decay can be expressed in terms of kinematic vari- 
ables and multipole matrix elements, using the formal- 
ism in [23] . This formalism is quite general, as it includes 
all multipoles and leaves the local weak current unspec- 
ified - so that exchange-current or other corrections to 
the usual one-body operators can be included. The tran- 
sition of interest involves an angular momentum change 
from 3=2 to J=0, with no change in parity. The relevant 
operators, expanded in the long- wavelength limit (since 
qR << he), are given in terms of the vector and axial 
vector currents, J(x) and J5(x) [23|. These operators, 
M2{q), T2^{q), and T^^^^{q), are the rank- two charge and 
transverse electric projections of the vector current and 
the transverse magnetic projection of the axial current, 
respectively. 

The vector current appearing in T^' ^i^s a one-body 
part and two-body exchange-current corrections, both of 
which are of the order (wnuc/c). A direct calculation of 
the important exchange current corrections would be dif- 
ficult to do precisely in a model calculation. Fortunately, 
under CVC and in the long wavelength limit, it is possi- 
ble to rewrite matrix elements of T2 ' in terms of those of 
M2 through Siegert's Theorem [2^. The benefit of this 
substitution is that the one-body charge operator is of or- 
der (wnuc/c)", while two-body corrections to this operator 



are of order (wnuc/c)^. Thus Siegert's theorem allows one 
to deal with a simpler operator that can be well approx- 
imated by its one-body form. The beta-decay rate then 
can be written in terms of M2{q) and T^'^^^{q). As dis- 
cussed below, the matrix element of the Coulomb opera- 
tor can be taken from the isobaric analogue gamma decay 
transition. Thus the nuclear-structure dependence of our 
result is effectively isolated in a single matrix-element ra- 

tio, (fr^^)/(M2). 

B. Calculation of the density matrix 

For weak and electromagnetic interactions, the inter- 
action can be well-approximated by considering only the 
one-body interaction [25|. We define |a) and \ j3) as com- 
plete single-particle wavef unctions, and we also define 
||a|) and as the reduced (in angular momentum and 
isospin) forms of these single-particle wavefunctions, i.e. 
jn, Z, J, t). We can then write a many-body multipole re- 
duced matrix element {J fTf\\(D jT\\JiTi) as the sum of 
single particle matrix elements (laHjOjTll multiplied 
by the density matrix elements ipafi, as defined in ^25j 

For the ^B(2+) ^Be(0+) transition, the coefficients 
for the density matrix were obtained from p-shell wave- 
functions calculated using the GLASGOW shell- model 12711 
with the Cohen-Kurath (8-16)2BME interaction [2J. 
The numerical values are given in Table [TTl 

TABLE II: Calculated density matrix elements for Ip-shell 
transitions for *B(2"*') *Be(0''') decay, using the Cohen- 
Kurath (8-16)2BME interaction [2^ and calculated using 
code obtained from [2^. The notation a ^ h denotes the 
density matrix element ipab, which connects the single parti- 
cle states with ji = a and jf = b. 



transition 


density matrix coefficient 




-0.206564 


±0.168713 


2 ~* 2 


-0.608563 



C. Single-particle matrix elements 

The single-particle operators relevant to this calcula- 
tion are defined in terms of spherical Bessel functions 
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jj and scalar and vector spherical harmonics Yj 



Mj 



and isoscalar components gives two results: either i?(i?2) 



Mf'(gx)=jj(gx)rf-'(17,) 



Mj 



i:f''{q^)=jj{qx)y%{n^)-CT. 



(1) 
(2) 



The matrix elements of these operators can be ex- 
panded in a harmonic oscillator basis and expressed as 
functions of the form I25l: 



(3) 

where Tj{qx.) is a generic single-particle operator reduced 
in spin/isospin, y = {^qb)^ , b=1.787 fm is the oscillator 
parameter [30|, and p(y ) is a polynomial of finite order 
in y, as tabulated in ■ 

Folding the matrix elements of [l^ with the Cohen- 
Kurath density matrix yields 



{Jf\\M2U] 



-0.23Fi 



magS I 



J,) = 



'47r 
0.04^:4^6 



^ ) exp(-2/), (4a) 



47r 



(-1) 



exp(-y), (4b) 



where Fi — 1.0 and Fa — 1-26 are the single-nucleon 
charge and axial- vector couplings. 

While the M2 and T^^s5 -jjiatiYx elements have iden- 
tical momentum dependence, a suppression of the T™^s5 
matrix element is caused by the cancelation to 20% in 
the sum of the density matrix elements: 



+ V'is = 0.17-0.21 



-0.04. 



(5) 



The uncertainty of this model dependence is difficult to 
quantify, and therefore we unfortunately cannot be as- 
sured that this result reflects the true strength of the axial 
contribution to this decay. Despite this, we are now able 
to obtain the beta decay rate for the ^B(2+) -> ®Be(0+) 
decay. We find that the vector contribution dominates, 
and the resulting branching ratio is 6.1 x 10~^. 



IV = 

O.OOiO.OSe^fm^ iindB{E2)is = 0.14±0.03e^fm4, or al- 
ternatively, B{E2)iv = 0.14±0.03e2fm'' and B{E2)is = 
0.00 ± 0.03 fm''. We will henceforth use the maximum 
(0.14 + 0.03 = 0.17e^fm"') and minimum (0.0) values for 
B{E2)\Y to constrain the beta decay transition rate. 

The rate of the nuclear electromagnetic transition rate 
in terms of the multipole operators defined in the previ- 
ous section is uM. 



2J. + 1 



(6) 



where lOj is the transition rate per unit time, kj is the 
gamma-ray momentum, a is the fine structure constant, 
J is the angular momentum change in the transition, 
and Ji is the initial spin of the nucleus. As with the 
beta decay expression, we can rewrite this in terms of 
the M2 matrix element using Siegert's Theorem [2^ and 
substituting J — Ji = 2, 



3 Sira {E, 
2~5 



Ef)' 



(7) 



The energy released in the isovector decay is 16.80 McV 
[2^ , and the corresponding decay rate from Eq. [7] yields 
a B{E2) value of 0.55 fm* for the transition, using 
the M2 matrix element obtained in the previous sec- 
tion. In order for this calculation to match the mea- 
sured B{E2) — 0.17 fm'*, the matrix element must be 
reduced by a factor of 1.8. To incorporate the second 
{B{E2)\Y ~ 0) constraint, we simply set M2 = in the 
beta decay calculation. 

If we assume the maximum possible strength of the 
vector matrix element from the gamma decay measure- 
ment, we obtain a transition rate of 1.8 x 10~^ s"'^. On 
the other hand, if we assume that the vector matrix ele- 
ment does not contribute, the transition rate is reduced 
to 9.1 X 10^^ s^'^. The respective branching ratios for 
these cases are 2.0 x 10~^ and 1.0 x 10^^. 



CONCLUSIONS 



D. Vector contribution derived from 
electromagnetic rate 

Since our shell model calculation showed dominance 
of the vector contribution to the decay, it is useful to 
compare that result with the value implied by measure- 
ments of isobaric analogue electromagnetic transitions. 
The doublet in ^Be at 16.6/16.9 MeV is composed of a 
T = state mixed with the T — \ isobaric analogue 
state of the *B ground state. The reduced strengths 
for the E2 transition from these states to the ground 
state of ^Be are B{E2)iqs = 0.068 ± 0.024 e^fm'' and 
B(£'2)i6.9 = 0.075 ± 0.013 e^fm'i [H. Due to ambi- 
guity in the phase of the matrix elements, decompos- 
ing the measured transition strengths into isovector and 



An accurate measurement of or limit on the *B(2+) 
^Be(0"'") beta decay branching ratio is of importance 
for the current generation of solar neutrino measure- 
ments. To the best of our knowledge, this work is 
the first published attempt at the measurement of the 
^B(2+) — > ^Be(0+) decay branching ratio. 

Our measurement of the ^B(2+) ^Be(0''") beta decay 
transition centers on the detection of the two low-energy 
alpha particles from the breakup of *Be. A significant 
challenge to this measurement is the characterization of 
the response of our silicon detector to low-energy alpha 
particles. In light of this, we have done two experiments 
to determine our detector response. The first measure- 
ment involved directly implanting alpha particles with 
kinetic energy between 86 and 453 keV into our detector. 
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while simultaneously measuring the detector deadlayer. 
The second measurement involved the implantation of 
-'^^N, which is unstable to beta-delayed 3a decay. The 
results from both of these measurements verify that the 
tabulated nuclear stopping powers for alpha particles in 
silicon are reasonably valid for the range of energies we 
are considering. However, the directly-implanted low- 
energy alpha particle measurements indicates that the 
uncertainty on the tabulated values may be larger than 
predicted. 

Using the measurements and techniques described 
above, we have measured the ^B(2+) *Be(0+) decay 
branching ratio to be (2.6 ± 3.3) x 10"^. This measure- 
ment is consistent with zero and can be used to place 
an upper limit on the branching ratio of 7.3 x 10~^ at 
the 90% confidence level. At this level, the ^6(2+) 
®Be(0"'") decay branch is not a significant background to 
measurements of the spectral shape of solar neutri- 
nos and the solar hep neutrino flux. Our calculation of 
the ®B(2+) ®Be(0"'") decay rate gives a branching ra- 
tio of 6.1 X 10~^, using values for the vector and axial- 
vector matrix elements obtained from the shell model. 
The model-dependent axial term in our calculation may 
be quite uncertain. We find a branching ratio range of 



1.0 X 10~^ to 2.0 X 10~^, based on the values for the vector 
matrix element obtained from the experiment of Ref. 
along with the axial- vector matrix element obtained from 
the shell model. The values are all consistent with our 
measured value. 
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